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ABSTRACT 

We report on the determination of the astrometric, spin, and orbital param¬ 
eters for PSR J1953+1846A, a “black widow” binary millisecond pulsar in the 
globular cluster M71. By using the accurate position and orbital parameters 
obtained from radio timing, we identified the optical companion in ACS/Hubble 
Space Telescope images. It turns out to be a faint (rripeoew > 24, nipsuw 23) 
and variable star located at only ~ 0.06" from the pulsar timing position. The 
light curve shows a maximum at the pulsar inferior conjunction and a minimum 
at the pulsar superior conjunction, thus confirming the association with the sys¬ 
tem. The shape of the optical modulation suggests that the companion star is 
heated, likely by the pulsar wind. The comparison with the X-ray light curve 
possibly suggests the presence of an intra-binary shock due to the interaction 
between the pulsar wind and the material released by the companion. This is 
the second identification (after COM-M5C) of an optical companion to a black 
widow pulsar in a globular cluster. Interestingly, the two companions show a 
similar light curve and share the same position in the color magnitude diagram. 
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Subject headings: Pulsars: Individual: PSR J1953+1846A, Globular clusters: 
Individual: M71 (NGC 6838), Techniques: photometric 


INTRODUCTION 


Millisecond pulsars (MSPs) are rapidly spinni ng neutron stars (NSs), formed in a binary _ 

syste m where, according to the canonical scenario flAlpar et al.lll982l:lBhattacharva &; van den Heuvel 
1 991 1), a slowly rotating NS is spun up through mass accretion from an evolving companion 
star. The accretion phenomena are usually observed during the phases of low mass X-ray 
binaries (LMXBs), which can be considered as the progenitors of MSPs. The process usu¬ 
ally leads to a deep transformation of both the accreting and the donor stars: the former 
is accelerated to millisecond spin periods, while the latter can evolve into an intermediate 
anomalous evolutionary phase (see, e.g., MSP-A in NGC 6397; Ferraro et ah 2001 al). before 
reaching the final stage of a (possibly He) white dwarf (WD; e.g. MSP-A in NGC 6752; 

Ferraro et al. 2003a). 


Although the Galaxy is ~ 100 times more massive than the entire Galactic GC ecosys¬ 
tem, about 40% of the known MSP population is found in GCs. Such an over-abundance is 
indicative of a strongly enhanced dynamical activity in these dense stellar systems. In fact, 
in the Galactic field the most plausible channel for the formation of MSPs is the evolution 
of primordial binaries, while in GCs dynamical interactions promote the formation of a con¬ 
spicuous number of exotic objects, such as blue straggler stars, X-ray binaries, cataclysmic 


variables and MSPs ( Bailvru Il992l: ICool et al.l Il995l : iFerraro et al.lll995l . l2001bl: IPoolev et al. 


20031: iRansom et al.l l2005aT) , which can be used to trace the complex i nterp l ay between dy 


namics and stellar evolution (e.g. 

Gooc 

man & Hut 

1989; 

Hut et al. 

1992 

Ferraro et al. 

1995; 

Possenti et al. 

2003; 

Ferraro et al. 

2003b, 

2009, 

2012). 


in the very centre of these systems, we expect to find a large number of NSs which are (or have 
been) affecte d by dyn amical processes such as tidal captures or exchange interactions (see 


e.g. Ivanova et al. 2008). In this respect, the study of optical companions to binary MSPs 
in GCs is of outmost importance since it opens the possibility to evaluate the frequency 
and time-scales of dynamical interactions in dense stellar systems, to explore the impact 
of dynamics on MSP evolutionary paths and to investigate stellar evolution under extreme 


1 Based on observations collected with the NASA/ESA HST (Prop. 12932), obtained at the Space Tele¬ 
scope Science Institute, which is operated by AURA, Inc., under NASA contract NAS5-26555. 
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conditions (see e.g. 


Ferraro et al. 

2003c; 

Sabbi et al. 

2003a 

b; 


Mucciarelli et ah 2013) 


Although the majority of binary MSPs have low-mass He WD companions, recent PSR 
searches have considerably increased the number of non-canonical binary MSPs. Among 
these, “black widows” (BWs) and “redbacks” (RBs) are of particular interest due to the 
presence of radio eclipses of the MSP signal, caused by a significant amount of ionized ma- 
terial ablated from the companion star because of the energy injected by the pulsar (PSR; 
e.g. Ruderman et ah 1989 ). The eclipsing regions are usually larger than the companion 


confirmed by several optical identifications (e.g. 

Edmonds et al. 

2002; 

Ferraro et al. 

2001a; 

Cocozza et al. 

2008; 

Pallanca et al. 

2010. 

2013, 

2014a 

). These systems are characterized 


by small orbital eccentricities, tight orbits (orbital periods P or b < 1 day) and small mass 
functions, thus implying the presence of a low mass companion. Indeed, RB companion 
stars have usually masses of ~ 0.1 — 0.5M 0 , while BW companions are much less massive 
(M < O.1M 0 ). Such a small value could be due to vaporization from the strong MSP ra¬ 
diation and relativistic wind; thus these systems may provide a possible explanation to the 
existence of isolated MSPs, even if their observed number is too large if com pared to the 
expected timescale for the total ablation of the companion stars (lEichler fe Levinson] 1988). 
The physical mechanism favoring the formation of a RB instead of a BW is not well un¬ 
derstood yet. Benvenuto et ah ( 120141 ) argued that RBs may evolve into BW systems, as a 
direct consequence of ablation processes, although this cannot apply to all RBs since some 
of them could evolve into canonical He WD systems. On the other hand, Chen et al. (2013) 
suggested that BWs cannot result from the evolution of RBs and that the discriminant factor 
leading to the formation of a RB instead of a BW is the reprocessing efficiency of the MSP 
emission by the companion star, likely related to geometrical factors. 


Since the first eclipsing MSPs were preferentially found in G Cs, i t was c ommonly be¬ 
lieved that these systems form exclusively in crowded environments (King et a 02QQ3 h mostly 
as a consequence of exchange interaction, and the few objects discovered in the Galactic field 
were born in GCs and later ejected. Indeed, at that time, the ratio between the number 
of eclipsing MSPs and canonical MSPs was significantly higher in GCs than in the Galac¬ 
tic field. This scenario has been recently altered by the di scovery of a large number of 


eclipsing MSPs in the G alactic field, both in blind surveys (IBureav et al.l 12006 


Bates et al. 


201 ll; Keith et al.ll2012l). and espe c ially in surveys ai med at the radio identification of Fermi 


sources 


e.g. Ransom et al. 2 0111: Ke ith et al . 201 ll). Hence, these objects can also form 


in the field from the undisturbed evolution of LMXB systems, with no need for dynami¬ 
cal interactions with other stars. In such a scenario, the fraction of eclipsing to canonical 




























































































4 


MSPs should be similar in the held and in GCs, irrespectively from their in teracti on r ate 
per binary, since binaries in very tight orbits are unlikely to be disrupted (Verbunt & Freire 


2 0141 1. As an example, the BWs J1518+0204C and J1807—2459A are located respectively 


in a GC with a very low interaction rate (M5) and in the one with the largest rate (NGC 
6544; iPallanca et al.l l2014ar I Verbunt fe Freirel 120141 1 . Still, both the objects appear to be in 
a relatively undisturbed systems (no orbital eccentricity). Understanding the formation of 
BWs and RBs, including the possible role of stellar interaction, strongly motivates multi- 
wavelength studies, both in GCs and in the Galactic held. 


Despite the importance of MSP optical companions, finding them in crowded stellar 
systems like GCs is extremely challenging. Only nine companions have been discoyered so 


Sia'urdsson el 

ah 2003), as expected from the cano 

— — — — — — — -l — -17 

meal formation scena 

i- 1 i— — n 

rio, five are RBs com- 

panions (see 

Ferraro et ah 2001a; Edmonds et ah 

2002; Cocozza et ah 

20 

08; Pallanca et ah 

2010. 20131 and onlv one is a BW companion (see Pallanca et ah 2014a 

). Here we report 


the radio timing ephemeris of PSR J1953+1846A (hereafter M71A) in the GC M71 (NGC 
6838) and the identification of its companion star. 


M71A is the only MSP known so fa r in M7 1 ( Ransom et al.ll2005bl) . M 71 is a low density 
GC ( log Pn = 2.83 in units of L 0 pc -3 ; lHarrislll996l) . in a disk-like orbit (jGeffert &; Maintz 


2000), loc ated at ~ 4 kpc from the Earth. It is one of the most metal-rich clusters among 


halo G Cs ( Harris 1996 ) and its surface brightness prohle shows an extended core (r c = 37.8"; 


Harr is 19961 ) and no signatures of core collapse. M71A was discovered in a ta rgeted survey 


of all GCs visible with the 305-m Arecibo radio telescope (IHessels et ah 120071 ). It is located 
at a = 19 h 53 m 46.42 s ; 8 = 18°47 / 04.84", at a pr ojected distance of only 20" (0.53 core radii) 
from the cluster center (iGoldsburv et all 201(1 ). and it has a spin period of ~ 4.9 ms and 
a low eccentricity orbit of ~ 4.2 hours. M71A is classihed as a BW. In fact, because of its 
very low mass function (/ = 1.6 ■ 1O _5 M 0 ), the companion is expected to have a minimum 
mass of ~ O.O32Af 0 . Moreover, as commonly found for BW systems, the radio signal shows 
eclipses for about 20% of the orbital period (at 1400 MHz observing frequency), likely due 
to stripped material from an evaporating companion. A Chandra X-ray observation of this 
cluster revealed a source in a position compatible with the PSR location arid a luminosity 
of about 10 31 ergs s _1 in the 0.3 — 8.0 keV spectral range ( Eisner et ah 2008h . The light 
curve is consistent with a non-steady source and the photon index (T = 1.89 ±0.32) suggests 
magnetospheric radiation and/or an emission from intra-binary shocks. In the context of an 
optical study of the M71 X-ray sources, IHuang et ah (J2010) suggested as possible optical 
companion to M71A a star located at ~ 0.1" from the radio PSR and lying along the red 
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side of the cluster main sequence (MS), in a region commonly occupied by binary systems. 
Nonetheless, its absolute magnitude {My ~ 8.5) implies a mass of about O.5M 0 , inconsistent 
with radio-derived mass function (in fact such a large mass would be compatibleyonly with 
a nearly face-on orbit, where no radio eclipses are expected). Hence, Huang et al . (120101) 
concluded that this object is unlikely to be the real companion, which could be still below 
the detection threshold or, alternatively, that M71A could be a hierarchical triple system. 


In Section 2 the radio timing analysis is presented, while Section 3 is devoted to the 
optical identification of the companion star. In Section 4 we discuss the results and in Section 
5 we summarize the current knowledge about this BW system. 


2. RADIO TIMING 


Timing observations were carried out with the 305-rn William E. Gordon telescope at the 
Arecibo Observatory in Puerto Rico, between MJDs 52420 (2002 May 26) and 53542 (2005 
June 21), with the initial discovery observations on MJD 52082 (2001 June 22) incorporated 
into the timing solution. The Gregorian L-band Wide receiver was used for the observa- 
tio ns, se nding dual-polarization data to the Wideband Arecibo Pulsar Processor (WAPP; 


see 


Dowd et ah 20001) autocorrelation spectrometers. For most of the timing observations, 3 


WAPPs were used, at frequencies centered near 1170, 1420 and 1520 MHz, although in the 
beginning, and occasionally thereafter, only one WAPP was used, and sometimes 4 WAPPs 
were used. The WAPPs were configured to autocorrelate 3-level samples with 256 lags and 
accumulate these for 64 /rs, then sum polarizations and write the results to disk as 16-bit 
numbers. See Hessels et al.l fl2007h for details of the observations. 

Offline, the PRESTO software^ was used to partially dedisperse the data into 32 sub¬ 
bands, to reduce the data size while still facilitating searches for further PSRs in the cluster 
(no additional radio PSRs were found, despite careful searches of the majority of the acquired 
timing data). The subbands were then folded modulo the best-known PSR ephemeris. A 
Gaussian profile was fit to the summe d p rofile from several observations for use as a standard 
profile, and the FFTFIT algorithm (ITavlorlll9921) was used to determine pulse times of arrival 
(TOAs). Time segments corresponding to eclipses and to the times when M71 transited at 
Arecibo (during which the telescope could not track the cluster) were not considered in the 
timing analysis. 


2 github.com / scottransom/presto 

















Timing analysis was performed with the TEMPO software packag<{] using the DE421 
Solar System ephemeris and the TT(BIPM) clock standard. The BT timing model of 
Blandford fe Teukolskvl (119761 ) was used, as the orbit has no significant eccentricity. The 
timing parameters are listed in Table |Tj and residuals are presented in Figure [H The root- 
mean-square postfit residual is 35 /rs. The reduced-y 2 of the fit is 4.8; however, we list the 
parameter uncertainties as reported by TEMPO without scaling, as the epoch-to-epoch 
wander in the residuals is likely due to the interactions between the two stars (see Figure [2]) 
rather than to any misestimation of the TOA uncertainties. The high-precision radio timing 
position is slightly offset (0.06") from the position of the optical counterpart (see Section [3]), 
bu t agrees w ithin the much larger uncertainty (0.2") of the latter. Following the reasoning in 
Freire et al.l (2005), we find the maximum possible acceleration due to the gravitational field 
of the cluster for this line of sight to be ±3.2 x 10~ 10 ms~ 2 . This implies that most of the 
observed pulse period derivative (P) is intrinsic. Further corrections due to the differential 
acceleration in the Galaxy (e.g. iNice fe Tavlorlll995t iReid et al.l 120141) are sm all. Given the 
small velocity dispersion in the core of the cluster (2.3 km s~ x ; Harris 199 6). the velocity 
of the PSR relative to that of the cluster should be very small; therefore its proper motion 
should be very similar to the proper motion of the cluster as a whole. A recen t measurement 
of the proper motion of the cluster amount s to 3. 93 mas yr -1 ( Kharchenko et al.l 20131), 
making the corresponding correction to P (IShklovskiil 19701) about half the size of that due 
to the Galactic acceleration. The timing data do not allow us to derive a reliable proper 
motion for the PSR. We use the range of allowed accelerations to constrain the intrinsic P 
as well as the characteristic age and surface magnetic field in Table [0 

The PSR is asymmetrically eclipsed between approximate orbital phases of 0.18 and 
0.35, where orbital phase 0.25 represents superior conjunction. The eclipses begin fairly 
abruptly but when the signal returns, it at first suffers excess dispersive delay due to ionized 
material within the orbit (Figured]). A discussion of the companion size and inclination angle 
are presented in Section 4. The mass loss from the companion star has a further manifestation 
in the variation of orbital parameters: Figure [2] shows the value of the time of ascending 


node passage for overlapping subsets of the data. The variation is comparable to that seen 

in other black widow eclipsing systems (e.g., 

Arzoumanian et ah 

1994 

: Ng et ah 

202) and 

significantly less than what is typically present in the redback systems 

(e-g-, 

Archibald et ah 


2 013 1). which have much more massive, likely non-degenerate companion stars. 


3 tempo. sourceforge. net 
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3. OPTICAL PHOTOMETRY OF THE COMPANION STAR 
3.1. Observations and data analysis 

The identification of the companion to M71A has been performed through two datasets 
of high resolution images acquired with the Wide Field Camera (WFC) of the Advanced 
Camera for Surveys (ACS) mounted on the Hubble Space Telescope (HST). The primary 
dataset has been obtained on 2013 August 20 (G012932, P.I.: Ferraro) and consists of a set 
of ten images in the F606W filter (with exposure times: 2 x 459 s; 3 x 466 s; 5 x 500 s) and 
nine images in the F814W filter (5 x 337 s; 3 x 357 s; 1 x 440 s). We also analyzed an archival 
dataset, obtained on 2006 July 1 (G01775, P.I.: Sarajedini) with the same instrument and 
same filters. It consists of four F606W images with an exposure time of 75 s and four F814W 
images with an exposure time of 80 s. 


The standard photometric analysis (see Dalessandro et al. 2008a,b) has been performed 
on the “flc” images, which are corrected for flat held, bias, dark counts and charge trans¬ 
fer efficiency. These images have been further corrected for “Pixel-Area-Map ”0 with stan¬ 
dard IRAF procedures. By using the DAOPHOT II ALLSTAR and ALLFRAME packages 
( Stetson 19871 ). we performed an accurate photometric analysis of each image. First of all, 
we modeled the Point Spread Function (PSF) by using a sample of ~ 200 bright but not 
saturated stars. The model has been chosen on the basis of a y 2 test and, in every image, 
the best fit is provided by a Moffat function flMoffatl 11969 ). Then we performed a source 
detection analysis, setting a 3a detection limit, where a is the standard deviation of the 
measured background. Once a list of stars was obtained, we performed a PSF-htting in each 
image. In the resulting catalog we included only stars present at least in half the images for 
each Liter. For each star, we homogenized the magnitudes estimated in different images and 
their weighted mean and standard deviation have been fina lly ad opted as the star mean mag¬ 


nitude and its related photometric error (see [Ferraro et ah 199 1. 199 2). However, in order to 


perform variability studies, for each source we also kept the homogenized magnitude mea¬ 
sured in each frame in both Liters. Then, instrumental magnitudes have been calibr ated t o 
the VEGAMAG system cross-correlatingp our catalog with that by [Anderson et ah ( 20081 ), 
using the ~ 7600 stars in common. 


4 For more details see the ACS Data Handbook. 

5 We used CataXcorr, a code which is specifically developed to perform accurate astrometric solutions. 
It has been developed by P. Montegriffo at INAF- Osservatorio Astronomico di Bologna. This package is 
available at http://davide2.bo.astro.it/^paolo/Main/CataPack.htmf and has been successfully used in a 
large number of papers by our group in the past 10 years. 






















Since the WFC images suffer heavily from geometric distortion, we co rrected the instru¬ 
mental positions ( x,y) by app lyi ng the equations reported by iMeurer et al. (120031) and using 
the coefficients in Hack & Cox (120011 ). Then we transformed instrumental posit i on s into 


th e ab solute astrometric system (ct, 5) using the stars in common with the Anderson et al. 
( 2008! ) catalog. The resulting astrometric solution has an accuracy of ~ 0.14" in a and of 
~ 0.13" in 5, corresponding to a total position accuracy of ~ 0.2". 


3.2. The companion to M71A 


The search for the companion star to M71A was performed by means of an accurate 
photometric analysis of all the detectable objects within a 5" x 5" wide region centered on 
the nominal position of the MSP. Figure [I] shows the zoomed (0.5" x 0.5") central part of 
that region. As can be seen, a relatively bright object is found to have a position compatible 
with the X-r ay sour ce (dashed circle) and the radio source (solid-line circle). This is the star 
proposed by Huang et al. ( 2Qlol ) to be the optical counterpart to M71A. However, a much 
fainter object, showing a strong variability, is visible in the figure. This is a quite promising 
object and it is located at a = 19 h 53 m 46.4062 s ; 6 = 18°47 , 04.793", only 0.06" from the 
radio position and 0.13" from the X-ray source, thus in perfect positional coincidence within 
our positional uncertainty (~ 0.2"). In the primary dataset, it has been detected in 9 (out 
of 10) images in the F606W filter, with a magnitude variation ranging from mp 606 W ~ 24.3 
to Hip 606 W ~ 27, while in the F814W filter has been detected in 6 images (out of 9) and the 
magnitude varies from m F 8 i 4 W ~ 23.4 to m F 8 i 4 W ~ 24.9. Unfortunately, the images in the 
archival dataset are too shallow to properly detect this faint object: in fact it turned out 
to be above the detection threshold in only one exposure in the F814W filter. For the four 
deep exposures of the primary dataset in which the star is not visible, we estimated an upper 
magnitude limit by simulating an artificial star of decreasing magnitude at the position of 
the candidate companion. The derived detection threshold turned out to be m F60 6w ~ 26.5 
and m F 8i4w 25.9. 


In order to reliably establish that the detected star is the binary companion to M71A, 
we built the light curve in both the available filters by folding the optical measurements 
with the orbital period and the ascending node time of the PSR (see Table [1]). The results 
are shown in Figure [Sj and in Table 2 we report the MJD of the images with their related 
orbital phases and magnitudes. As can be seen, the light curves show a sinusoidal modula¬ 
tion spanning at least three magnitudes and it is fully consistent with the orbital period of 
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the binary system. This establishes the physical connection between the variable star and 
the MSP. Indeed the exposures in which the star is not detected nicely correspond to the 
light curve minima. The curves have a maximum at <f> ~ 0.75, corresponding to the PSR 
inferior conjunction (where we observe the companion side facing the PSR) and a minimum 
at (j) « 0.25, corresponding to the PSR superior conjunction (where we observe the back side 
of the companion). This behavior is indicative of a strong heating of the companion side 


ties of other similar objects (e.g. 

Stapoers et al 

. 2001; 

levnolds et al. 

2007; 

Pallanca et al. 

2012 ; 

Breton et al. 

2013; 

Pallanca et al. 

2014a; 

Li et al 

2014 

). For the sake of comparison. 


in Figure [6] we plot the light curve (folded following the same procedure described above) of 
the possible companion suggested by Huang et al. fl2010l ). As can be seen the star does not 
show any significant flux variation. 


All these pieces of evidence suggest that the faint variable (which we name COM-M71A) 
is the optical companion to the BW M71A. It is the tenth MSP optical companion and the 
second to a BW system in a GC. In the color-magnitude diagram (CMD), COM-M71A is 
located at faint magnitudes in a region between the MS and the WD cooling sequences, 
where no normal GC stars are expected. This position is indicative of a non degenerate or 
semi-degenerate, low-mass and swollen star. Interestingly, the position of this object in the 
CMD is quite similar to that of COM-M5C, the onl y companion to a BW system known in 
GCs up to now and recently identified by Pallanca et. al, (2014a) in the GC M5. 


4. DISCUSSION 

Since the available data do not uniformly sample the orbital phases of the system in 
either the F606W or the F814W filters (see upper and middle panels of Figure [5]), in order 
to accurately determine the light curve of the companion star we combined the two datasets 
together, by applying a 0.95 mag shift to the F814W magnitudes (bottom panel of Figure [5]). 
We then used the software GRATIS^ and a x 2 criterion to determine the two harmonics best 
fit mode|j]to the curve. This is shown as a solid line in the bottom panel of Figure [5] Finally, 


6 “Graphical Analyzer for Time Series” is a software aimed at studying stellar variability phenomena. 
Developed by Paolo Montegriffo at INAF-Osservatorio Astronomico di Bologna. 

7 Note that a single harmonic model (i.e. a sinusoidal function) does not provide a good match of the 
observed light curve. 
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we verified that the same solution also provides a good fit to the light curves in each filter 
separately. Indeed, the reduced y 2 turned out to be 1.50 for the F606W filter, and 1.75 in 
F814W (see the solid curves in the upper and middle panels of Figure [5]). This demonstrates 
that no significant modulation of the stellar color (temperature) along the orbit is measur¬ 
able from the available dataset, and a much finer sampling of the light curve is needed to 
provide additional clues on this possibility. In Table 3 we report the maximum and minimum 
values for both the magnitude and the flux in each filter, evaluated from the best-fit model 
by following the procedure described in Bohlin (2012) for the ACS. The uncertainties are 
calculated by using the mean photometric errors for stars with similar magnitudes. The 
magnitude shift needed to superimpose the F814W light curve to that in the F606W filter 
implies a color index of 0.95 ± 0.12 for the companion star. By adopting a 0.54 Mp WD 
cooling sequence from the BaSTI catalog^ ( Manzato et al. 2008 : Salaris et al. 2010 1. this 


with those evaluated for other BW svstems (e.g. 

Stampers et al. 

2001; 

Pallanca et al. 

2012; 

Breton et al. 

2013; 

Pallanca et al. 

2014a; 

Li et al. 

2014 

)• 


In Figure [Tj we show the CMD, with the shaded rectangle marking the region occupied 
by COM-M71A during the orbital period. The height of the rectangle corresponds to the 
maximum niF 606 W magnitude difference expected from the best-fit model shown in Figure 0 
while the width corresponds to the photometric error at the minimum luminosity. As al¬ 
ready mentioned, the star is located between the MS and the WD cooling sequence, and it 
spans a range of about three magnitudes. Of particular interest is the predicted star position 
during the PSR superior conjunction, where we expect to see the stellar side not exposed 
to the PSR flux (clearly this is exactly the case only for i = 90°). The CMD position of 
COM-M71A in this phase could be compatible with the He WD cooling sequence, suggest¬ 
ing a semi-degenerate stellar structure. However, at these low luminosities, our analytical 
model is not appropriately constrained by data. Therefore, in order to confirm this possibil¬ 
ity, further observations are needed. In principle, the companion mass can be constrained 
from the comparison of its CMD position and theoretical isochrones. However, in the case 
of strongly perturbed stars the mass inferred in this way can be overestimated, or sugges- 
tive of inclination angles too small to be consistent with the pre sence of radio eclipses (see 


Ferraro et al.ll2003al: iPallanca et al.ll2010l; iMucciarelli et al.l 120131 ) . In our case, not only the 


companion position in the CMD is clearly out of the canonical evolutionary sequences, but 
also its minimum luminosity is not properly constrained by the observations. 


8 http: //basti.oa-teramo. inaf.it 
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Assuming that the companion optical emission is mostly due to black-body radiation, 
the luminosity and temperature of this star would be consistent with an object of radius 
Rbb < O.O2i? 0 . However, the companions to BWs usually suffer from strong tidal distortion 
due to the interaction with the PSR, therefore they are swollen up and possibly they can 
even fill their Roche Lobes. Furthermore, the presence of radio eclipses suggests that the 
simple Rbb is a gross underestimate of the true stellar radius. Indeed, the Roche_Lobe 


radius is far more appropriate to describe the size of the compan ion (e.g. IStappers et al. 


19961 : Pallanca et al.l 


2012 


(11983 1. the Roche Lobe radius can be computed as: 


Breton et al.ll2013t iPallanca et al.ll2014af) . According to [Eggleton 


R 


RL 


0A9qi 


0 .6<p + In (1 + qa 


( 1 ) 


where q is the ratio between the companion and the PSR masses and a is the orbital sepa¬ 
ration. Combining this relation with the PSR mass function, assuming a NS mass ranging 
from 1.2M 0 to 2.5 M & ( Ozel et al. 2012 ) and an inclination angle ranging from 0° to 90°, we 
find 0.22 R e < R rl < 1.24 R & . 


Interestingly, the light curve shape presents a hint of asymmetric structure in both Li¬ 
ters: the increase to the maximum seems to be smoother than the decrease to the minimum. 
Despite the low statistic, this behavior could be du e to a s light a synch ronous rotating com¬ 


panion, as in the case of PSR J2051—0827 (see Stappers et ah 2001). This could be the 


result of a tidal torque from the wind of a magnetically active star, which can result in a 
companion angular velocity that differs from the orbital angular velocity. Moreover, in this 
case the angular velocity could be subject to a variation with time due to a secular time 
dependence o f the orbital period, due itself to a variation of the companion quadrupole mo - 
ment (see e.g. lApplegate fe Shahamlll994l: iDoroshenko et al.ll200ll: iLorimer fe Kramerll20121 ) . 
However, in order to probe this intriguing possibility, an uniform sampling of the light curve 
from new observations is needed. 


4.1. Reprocessing efficiency and Roche Lobe filling factor 

Under the assumption that the optical magnitude modulation is mainly due to the heat¬ 
ing of the companion surface by the PSR flux, we can compare the observed flux amplitude 
of the light curve with the expected one (AF exp ) as a function of the inclination angle i , by 
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the following relation (jPallanca et al.ll2014al) : 


AF eX p ( Z ) T) 


E 


Rcom(^) 


eu 


dpsR 


( 2 ) 


where rj is the reprocessing efficiency under the assumption of a PSR isotropic emission, 
Rcom{ i) = f Rrl( i) is the companion star radius, where / is the volume-averaged Roche 
Lobe Filing fact or, dps n is the MSP distance, assumed to be equal to the GC distance 
(dpsR = 4.0 kpc; [Har ris 19961 ) and e(z) parametrizes the difference of the heated surface vis¬ 
ible to the observer between maximum and minimum, as a function of the inclination angle. 
E = 4-7Tis the PSR spin-down luminosity where / is the momentum of inertia. Using 
the the spin period and its intrinsic first derivative obtained from radio timing (Table 1) and 
assuming / = 10 45 g cm 2 , we found that 4.6 • 10 33 ergs s _1 < E < 5.8 ■ 10 33 ergs s _1 , typical 
values within the Galactic eclipsing MSP population. Setting AF exp = AF 0 b s in the F606W 
filter (see Table 3), we evaluated the reprocessing efficiency as a function of the inclination 
angle for different values of the Roche Lobe filling factor. Results are shown in Figure |8j 
As can be seen, for high inclination angles and a Roche Lobe filling companion, the repro¬ 
cessing efficiency is ~ 5%, while for Filing factor / = 0.8 is ~ 8%. A typical value of 15% 
(jBreton et al.l 120131 1 would be consistent with / ~ 0.6. Values of / < 1 would be plausible 
since some works showed that BW companions not always comple tely Fll their Roche Lobe 


(e-g- 

the F814W Flter. 


Callanan et al.lll995l: IStappers et al.lll999[ iBreton et al.l 120131) . Similar results hold for 


It is worth noting that by using R BB instead of R rl for the stellar radius, the efficiency 
increases over 100% for almost every meaningful conFguration. This can be admitted only 
if an anisotropic PSR emission is assumed. However, the presence of long radio eclipses and 
the behavior of similar objects is a strong indication that R BB heavily underestimates the 
stellar true radius. 


4.2. A comparison between M71A and M5C 


So far, the optical companion to PSR J1518+0204C (hereafter M5C) was the only BW 
companion known in a GC. M5C is a 2.48 ms PSR with an orbital period of ~ 2.1 hr located in 
the GC M5. Its radio timing a nd t he optical photometry of the companion star (COM-M5C) 
is discussed in Pa llan ca et ah (j2014a|). In section 3, we anticipated some interesting analogies 
between this star and COM-M71A. In order to further investigate similarities between these 
BW companions, we compared the optical properties of the two objects. Figure [9] shows 
their light curves, with the magnitudes reported to the absolute values. Indeed, despite 
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the low sampling of the COM-M5C light curves, these two objects seem to have a quite 
similar optical behavior. As reference, we used C0M-M71A analytical models (solid lines) 
for the COM-M5C, from which we inferred that a similar light curve structure could hold 
even for C0M-M5C, being more appropriate than the simple sinusoid (dashed lines) used 
by Pallanca et ah (l2014al) . given the sparse number of measurements that prevent them to 
build an accurate model. Figure flOl shows the position of the two objects in the CMD. Again, 
considering the uncertainties in COM-M5C magnitudes and colors, we found that they are 
located in the same region, suggesting a common evolutive path for these low-mass, possibly 
non degenerate, swollen and heated companions. Interestingly, in the CMD the two BW 
companion s are located in a region completely different from that usually occupied by RBs 
(jPallancal 2014b). Of course, additional identifications of BW companions are needed to 
firmly characterize the evolution of these objects. In addition, using equation j2J) for COM- 


M5C, s etting the filling factor / = 1 and using the spin-down period from Pallanca et ah 
(I2014al) to evaluate the spin-down luminosity (0.7 • 10 34 ergs s _1 < Emsc < 3 • 10 34 ergs s -1 ), 
we found a reprocessing efficiency rj ~ 5 — 20%, a value fully in agreement with what found 
for COM-M71A, thus further strengthening the analogies between these two systems. 


4.3. Comparing X-ray and optical light curves 


Usually, BWs with a high energy counterpa rt do not show any appreciable X-ray vari- 
ability related to their orbital period (see, e.g., iBogdanov et al.l 120061; iGentile et al.l 12014 
respectively for the BWs in the GC 47 Tucanae and in the Galactic held). However, this 
could be an observational bias, due to the lack of deep enough and systematic surveys of BWs 
in the X-rays. On the other hand, it is worth noting that several RB systems clearly show 


orbital X-ray modulation likely due to the presence of intra-binary shocks (IBogdanov et al . 


2006, 2011 


2 0141). M71A is an exception, since it has been found to show periodic X-ray 


variability (Eisner et ah 2008). Very interestingly, the determination of COM-M71A opti¬ 
cal light curve offers the opportunity to perform a comparison between the two. The most 
intriguing feature emerging from the comparison of the light curves (both folded with the 
binary system parameters) is that the phase spanned by the radio eclipses (0.18 < 0 < 0.35) 
does not correspond to the phase of the X-ray minimum (0 < 0 < 0.2), but it nicely lines 
up with the optical minimum (0 « 0.25). Thus we found that the X-ray minimum pre¬ 
cedes the optical PSR superior conjunction. A similar effect was already observed for the 


RB 47TucW, a 2.35 ms binary 

MSP with an orbital period of ~ 3.2 hr and a companion 

mass of ~ 0.15( 

Camilo et a 

2000 

). whose optical light curves indicate the presence of a 

strong heating (] 

Almonds et al. 

2002 ), 

as usually observed for BW systems. For this object 

Bogdanov et al. 

(2006) argue that the X-ray variability can be attributed to the presence 
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of an intra-binary shock that is eclipsed by the companion star. The length of the X-ray 
eclipse suggests that this shock is located closer to the companion star than to the MSP. 
In particular this behavior could be due to a swept-back shocked region, produced by the 
interaction between the PSR wind and the stream of gas issuing from the inner Lagrange 


point LI, elongated perpendicular to the semi-major axis of the binary (see Bogdanov et al. 


2005 . for a detailed description). Despite the low X-ray statistics, this is likely to be the 


case also for M71A, where the intra-binary shock could be eclipsed just before the PSR su¬ 
perior conjunction. Even for a companion that is under-filling its Roche Lobe, this shocked 
region can b e created t hanks to the stellar wind which can result in mass outflow through 
LI ( Bogdanov et al. 2005 ). 


As discussed in Bogdanov et al.J (2005J), the Accreting Millisecond X-ray Pulsar (AMXP) 
SAX J1808.4—3658, during quiescent states, shows several analogies with the RB 47TucW, in 
terms of both the X-ray spectrum and the optical variability. Based on the discussion above, 
these properties are also similar to those observed for M71A and, very interestingly, even 



these objects are urged to unveil connections between AMXPs and eclipsing MSPs, and 
between BWs and RBs. Indeed, several important new connections between AMXPs and 
RBs have been made in the last ye ars, esp ecially with the discoveries of system s transit ioning 
from one state to the other (see 


1 — 1 - —'.I 

Archibald et al. 

2009: 

Papitto et al. 

2013; 

Patruno et al. 


20141: iBassa et al.l 12014 : IStappers et al.ll2014j ). 


5. SUMMARY AND CONCLUSION 

We presented a phase-connected radio timing solution for the BW PSR J1953+1846A, 
which includes a very precise position and orbital parameters determination. Taking ad¬ 
vantage of this precise measure of the PSR position, we have used a set of high resolution 
ACS/HST images to search for the companion star in the optical bands. We identified a 
faint and strongly variable star (COM-M71A), showing a modulation of at least three mag¬ 
nitudes in both used Liters (F606W and F814W). In the CMD, COM-M71A lies in the 
region between the cluster MS and the WD cooling sequences, thus suggesting that it is 
a low-mass, non-degenerate or at least semi-degenerate star, with a temperature of about 
5100 K. Unfortunately, because of its faintness, it was detectable only in 16 out of 27 images, 
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mostly during the PSR inferior conjunction. The light curve shows a sinusoidal shape with 
a period fully consistent with the binary MSP. The maximum, during the PSR inferior con¬ 
junction, and the minimum, during the PSR superior conjunction suggest a strong heating 
of the companion star side exposed to the PSR flux. Such a behavior is in good agreement 
with that observed for similar objects in the Galactic held. By modeling the light curve, we 
showed that the companion reprocessing efficiency of the PSR energy is ~ 5% for a Roche 
Lobe filling companion, while a typical value of 15% is reasonable by assuming a filling fac¬ 
tor of 0.6. The comparison between the optical and X-ray light curves suggests the possible 
presence of intra-binary shocks, similarly to what observed for the RB 47TucW. A X-ray 
and optical follow-up will highlight the presence of this shocks and, possibly, will allow to 
characterize their properties and structure. Unfortunately, the star is too faint to allow a 
spectroscopic follow-up with the available instruments. However, an optimized photometric 
follow-up would provide the opportunity of better constrain the system properties, and by 
using, for example, phase-resolved observations with a narrow Ha filter we could constrain 
the presence of ionized material, eventually related to the intra-binary shocks. 

COM-M71A is, so far, the second BW optical companion identified in a GC after COM- 
M5C in M5. Interestingly, both the light curve shape and the position in the CMD are quite 
similar in the two systems. This suggests that probably the two objects undergo a similar 
evolutionary path. Even though the statistic is by far too limited to draw any meaningful 
conclusion, at the moment no significant difference with the BW optical companions ob¬ 
served in the held can be evidenced, probably suggesting that no dynamical interactions are 
strictly needed for forming these systems. 


6. Acknowledgement 

This research is part of the project Cosmic-Lab (http://www.cosmic-lab.eu) funded by 
the European Research Council under contract ERC-2010-AdG-267675. J.W.T.H. acknowl¬ 
edges funding from an NWO Vidi fellowship and ERG starting grant “DRAGNET” (337062). 


REFERENCES 

Alpar, M. A., Cheng, A. F., Ruderman, M. A., & Shaham, J. 1982, Nature, 300, 728 
Anderson, J., Sarajedini, A., Bedin, L. R., et al. 2008, AJ, 135, 2055 
Applegate, J. H., & Shaham, J. 1994, ApJ, 436, 312 

Archibald, A. M., Stairs, I. H., Ransom, S. M., et al. 2009, Science, 324, 1411 


16 


Archibald, A. M., Kaspi, V. M., Hessels, J. W. T., et al. 2013, arXiv:1311.5161 
Arzoumanian, Z., Fruchter, A. S., & Taylor, J. H. 1994, ApJ, 426, L85 
Bailyn, C. D. 1992, ApJ, 392, 519 

Bassa, C. G., Patruno, A., Hessels, J. W. T., et al. 2014, MNRAS, 441, 1825 

Bates, S. D., Bailes, M., Bhat, N. D. R., et al. 2011, MNRAS, 416, 2455 

Benvenuto, O. G., De Vito, M. A., & Horvath, J. E. 2014, ApJ, 786, LL7 

Bhattacharya, D., & van den Heuvel, E. P. J. 1991, Phys. Rep., 203, 1 

Bogdanov, S., Grindlay, J. E., & van den Berg, M. 2005, ApJ, 630, 1029 

Blandford, R., & Teukolsky, S. A. 1976, ApJ, 205, 580 

Bogdanov, S., Archibald, A. M., Hessels, J. W. T., et al. 2011, ApJ, 742, 97 

Bogdanov, S., Grindlay, J. E., Heinke, C. O., et al. 2006, ApJ, 646, 1104 

Bogdanov, S., Patruno, A., Archibald, A. M., et al. 2014, ApJ, 789, 40 

Bohlin, R. C. 2012, Instrument Science Report ACS 2012-01 

Burgay, M., Joshi, B. C., D’Amico, N., et al. 2006, MNRAS, 368, 283 

Breton, R. P., van Kerkwijk, M. H., Roberts, M. S. E., et al. 2013, ApJ, 769, 108 

Callanan, P. J., van Paradijs, J., & Rengelink, R. 1995, ApJ, 439, 928 

Camilo, F., Lorimer, D. R., Freire, P., Lyne, A. G., & Manchester, R. N. 2000, ApJ, 535, 
975 

Campana, S., D’Avanzo, P., Casares, J., et al. 2004, ApJ, 614, L49 

Chen, H.-L., Chen, X., Tauris, T. M., & Han, Z. 2013, ApJ, 775, 27 

Cocozza, G., Ferraro, F. R., Possenti, A., et al. 2008, ApJ, 679, L105 

Cool, A. M., Grindlay, J. E., Cohn, H. N., Lugger, P. M., & Slavin, S. D. 1995, ApJ, 439, 
695 

Dalessandro, E., Lanzoni, B., Ferraro, F. R., et al. 2008a, ApJ, 677, 1069 
Dalessandro, E., Lanzoni, B., Ferraro, F. R., et al. 2008b, ApJ, 681, 311 


-17- 

Doroshenko, O., Lohmer, O., Kramer, M., et al. 2001, A&A, 379, 579 

Dowd, A., Sisk, W., & Hagen, J. 2000, IAU Colloq. 177: Pulsar Astronomy - 2000 and 
Beyond, 202, 275 

Edmonds, P. D., Gilliland, R. L., Heinke, C. O., Grindlay, J. E., & Camilo, F. 2001, ApJ, 
557, L57 

Edmonds, P. D., Gilliland, R. L., Camilo, F., Heinke, C. O., & Grindlay, J. E. 2002, ApJ, 
579, 741 

Eggleton, P. P. 1983, ApJ, 268, 368 

Eisner, R. F., Heinke, C. O., Cohn, H. N., et al. 2008, ApJ, 687, 1019 
Eichlcr, D., & Levinson, A. 1988, ApJ, 335, L67 

Ferraro, F. R., Clementini, G., Fusi Pecci, F., & Buonanno, R. 1991, MNRAS, 252, 357 

Ferraro, F. R., Clementini, G., Fusi Pecci, F., Sortino, R., & Buonanno, R. 1992, MNRAS, 
256, 391 

Ferraro, F. R., Fusi Pecci, F., & Bcllazzini, M. 1995, A&A, 294, 80 

Ferraro, F. R., Possenti, A., D’Amico, N., & Sabbi, E. 2001a, ApJ, 561, L93 

Ferraro, F. R., D’Amico, N., Possenti, A., Mignani, R. P., & Paltrinicri, B. 2001b, ApJ, 561, 
337 

Ferraro, F. R., Possenti, A., Sabbi, E., & D’Amico, N. 2003a, ApJ, 596, L211 

Ferraro, F. R., Possenti, A., Sabbi, E., et al. 2003b, ApJ, 595, 179 

Ferraro, F. R., Sabbi, E., Gratton, R., et al. 2003c, ApJ, 584, L13 

Ferraro, F. R., Dalessandro, E., Mucciarelli, A., et al. 2009, Nature, 462, 483 

Ferraro, F. R., Lanzoni, B., Dalessandro, E., et al. 2012, Nature, 492, 393 
Freire, P. C. C., Hessels, J. W. T., Nice, D. J., et al. 2005, ApJ, 621, 959 
Fruchter, A. S., Stinebring, D. R., & Taylor, J. H. 1988, Nature, 333, 237 
Geffert, M., & Maintz, G. 2000, A&AS, 144, 227 

Gentile, P. A., Roberts, M. S. E., McLaughlin, M. A., et al. 2014, ApJ, 783, 69 



- 18 - 

Goldsbury, R., Richer, H. B., Anderson, J., Dotter, A., Sarajedini, A., & Woodley, K. 2010, 
AJ, 140, 1830, Erratum AJ, 142, 66 (2011) 

Goodman, J., & Hut, P. 1989, Nature, 339, 40 

Hack, W., & Cox, C. 2001, Instrument Science Report ACS 2001-008, 10 pages, 8 

Harris, W. E. 1996, AJ, 112, 1487 

Hessels, J. W. T., Ransom, S. M., Stairs, I. H., Kaspi, V. M., & Freire, P. C. C. 2007, ApJ, 
670, 363 

Huang, R. H. H., Becker, W., Edmonds, P. D., et al. 2010, A&A, 513, A16 

Hut, P., McMillan, S., Goodman, J., et al. 1992, PASP, 104, 981 

Ivanova, N., Heinke, C. O., Rasio, F. A., Bclczynski, K., & Fregeau, J. M. 2008, MNRAS, 
386, 553 

Keith, M. J., Johnston, S., Ray, P. S., et al. 2011, MNRAS, 414, 1292 

Keith, M. J., Johnston, S., Bailes, M., et al. 2012, MNRAS, 419, 1752 

King, A. R., Davies, M. B., & Beer, M. E. 2003, MNRAS, 345, 678 

Kharchenko, N. V., Piskunov, A. E., Schilbach, E., Roser, S., & Scholz, R.-D. 2013, A&A, 
558, AA53 

Li, M., Halpern, J. P., & Thorstensen, J. R. 2014, ApJ, 795, 115 

Lorimcr, D. R., & Kramer, M. 2012, Handbook of Pulsar Astronomy, by D. R. Lorimer , 
M. Kramer, Cambridge, UK: Cambridge LIniversity Press, 2012, 

Manzato, P., Pietrinfcrni, A., Gasparo, F., Taffoni, G., & Corclier, D. 2008, PASP, 120, 922 

Meurer, G. R., Lindler, D., Blakeslee, J. P., et al. 2003, HST Calibration Workshop : Hubble 
after the Installation of the ACS and the NICMOS Cooling System, 65 

Moffat, A. F. J. 1969, A&A, 3, 455 

Mucciarelli, A., Safaris, M., Lanzoni, B., et al. 2013, ApJ, 772, L27 

Ng, C., Bailes, M., Bates, S. D., et al. 2014, MNRAS, 439, 1865 

Nice, D. J., & Taylor, J. H. 1995, ApJ, 441, 429 



19 - 

Ozel, F., Psaltis, D., Narayan, R., & Santos Villarreal, A. 2012, ApJ, 757, 55 

Pallanca, C., Dalessandro, E., Ferraro, F. R., et al. 2010, ApJ, 725, 1165 

Pallanca, C., Mignani, R. P., Dalessandro, E., et al. 2012, ApJ, 755, 180 

Pallanca, C., Dalessandro, E., Ferraro, F. R., Lanzoni, B., & Beccari, G. 2013, ApJ, 773, 
122 

Pallanca, C., Ransom, S. M., Ferraro, F. R., et al. 2014, arXiv: 1409.1424 
Pallanca, C. 2014, arXiv: 1405.2898 

Papitto, A., Ferrigno, C., Bozzo, E., et al. 2013, Nature, 501, 517 

Patruno, A., Archibald, A. M., Hessels, J. W. T., et al. 2014, ApJ, 781, LL3 

Phinney, E. S. 1992, Royal Society of London Philosophical Transactions Series A, 341, 39 

Pooley, D., Lewin, W. H. G., Anderson, S. F., et al. 2003, ApJ, 591, L131 

Possenti, A., D’Amico, N., Manchester, R. N., et al. 2003, ApJ, 599, 475 

Ransom, S. M., Hessels, J. W. T., Stairs, I. H., et al. 2005a, Science, 307, 892 

Ransom, S., Hessels, J., Stairs, I., et al. 2005b, Binary Radio Pulsars, 328, 199 

Ransom, S. M., Ray, P. S., Camilo, F., et al. 2011, ApJ, 727, LL16 

Reid, M. J., Menten, K. M., Brunthaler, A., et al. 2014, ApJ, 783, 130 

Reynolds, M. T., Callanan, P. J., Fruchter, A. S., et al. 2007, MNRAS, 379, 1117 

Roberts, M. S. E., Mclaughlin, M. A., Gentile, P., et al. 2014, Astronomische Nachrichten, 
335, 313 

Ruderman, M., Shaham, J., Tavani, M., & Eichler, D. 1989, ApJ, 343, 292 
Sabbi, E., Gratton, R. G., Bragaglia, A., et al. 2003a, A&A, 412, 829 
Sabbi, E., Gratton, R., Ferraro, F. R., et al. 2003b, ApJ, 589, L41 

Safaris, M., Cassisi, S., Pietrinferni, A., Kowalski, P. M., & Isern, J. 2010, ApJ, 716, 1241 
Shklovskii, I. S. 1970, Sov. Astron., 13, 562 


-20 - 

Sigurdsson, S., Richer, H. B., Hansen, B. M., Stairs, I. H., & Thorsett, S. E. 2003, Science, 
301, 193 

Stappers, B. W., Bessell, M. S., & Bailes, M. 1996, ApJ, 473, L119 

Stappers, B. W., van Kerkwijk, M. H., Lane, B., & Kulkarni, S. R. 1999, ApJ, 510, L45 

Stappers, B. W., van Kerkwijk, M. H., Bell, J. F., & Kulkarni, S. R. 2001, ApJ, 548, L183 

Stappers, B. W., Archibald, A. M., Hessels, J. W. T., et al. 2014, ApJ, 790, 39 

Stetson, P. B. 1987, PASP, 99, 191 

Taylor, J. H. 1992, Phil. Trans. Roy. Soc. A, 341, 117 

Verbunt, F., & Freire, P. C. C. 2014, A&A, 561, AAll 


This preprint was prepared with the AAS IATpX macros v5.2. 



- 21 - 


Table 1. Timing parameters for PSR J1953+1847 (M7lA) a 


Parameter Value 


Measured Parameters 

Right ascension, a (J2000). 19 h 53 m 46^41966(3) 

Declination, S (J2000) . +18° 47' 04"8472(7) 

Spin frequency F, Hz . 204.57006473073(3) 

Spin frequency derivative, F (10 -15 ) (s -2 ) . —2.0299(3) 

Spin frequency second derivative, F (10~ 25 ) (s~ 3 ) 5.4(3) 

Epoch (MJD) . 52812.0 

Dispersion measure, DM (cm _3 pc) . 117.3941(15) 

DM derivative (cm _3 pcyr _1 ) . —0.0274(17) 

Orbital period, P b (d) . 0.1767950297(2) 

Projected semi-major axis, x (s) . 0.0782246(12) 

Epoch of Ascending Node, T 0 (MJD) . 52811.8761877(3) 


Derived Parameters 

Spin period, P (ms) . 4.8883007458412(6) 

Spin period derivative P (10 -20 ) . 4.8506(8) 

Spin period second derivative P (10 -29 ) (s _1 ) ... —1.28(7) 

Angular offset from cluster centre 9± (') . 0.33 

Intrinsic period derivative P int (10~ 2 °). 4.3 < P int < 5.4 

Characteristic age r c (Gyr) . 1.4 < r c < 1.8 

Surface magnetic field B 0 (10 8 ) (gauss) . 4.6 < B 0 < 5.2 

Mass function / (M 0 ) . 0.0000164427(8) 

Minimum companion mass m c (M 0 ) b . 0.032 


a Numbers in parentheses are uncertainties in the last digits quoted. 

b Computed assuming an orbital inclination angle of 90° and a pulsar 
mass of 1.4 M 0 . 




























Table 2: Optical observations of COM-M71A 


0 

t ( MJD ) 

m F606W 

111F814W 

0.02 

56524.57602385 

— 

24.9 ± 0.2 

0.06 

56524.58290459 

26.5 ± 0.2 

— 

0.10 

56524.59012681 

26.8 ± 0.3 

— 

0.22 

56524.43431532 

> 26.5 

— 

0.26 

56524.44193085 

26.8 ± 0.4 

— 

0.31 

56524.44960589 

— 

> 25.9 

0.33 

56524.63013255 

— 

> 25.9 

0.35 

56524.45750982 

27.1 ± 0.4 

— 

0.36 

56524.63586163 

— 

> 25.8 

0.40 

56524.64270200 

26.8 ± 0.2 

— 

0.44 

56524.64977385 

— 

24.87 ± 0.09 

0.48 

56524.65661404 

25.54 ± 0.08 

— 

0.50 

53867 . 78512088 “ 

— 

24.3 ± 0.2 

0.57 

56524.49569959 

— 

23.97 ± 0.09 

0.60 

56524.50166033 

— 

23.76 ± 0.09 

0.64 

56524.50885348 

24.44 ± 0.03 

— 

0.68 

56524.51627848 

— 

23.40 ± 0.05 

0.73 

56524.52347163 

24.30 ± 0.03 

— 

0.94 

56524.56203070 

— 

24.26 ± 0.06 

0.98 

56524.56891144 

25.6 ± 0.1 

— 


Note. — Orbital phases ((f>), corresponding MJD (t) of the observations and observed magnitudes or 
upper-limits in both filters. 

“ This is the only image of the archival dataset where the companion star is above the detection threshold. 
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Table 3: Optical properties of COM-M71A 



F606W 

F814W 

^bright 

24.31 ±0.01 

23.37 ±0.02 

® faint 

27.62 ±0.09 

26.7 ±0.1 

F bright (10 -17 ergs cm" 2 s" 1 ) 

126 ±1 

128 ±2 

F faint (10” 17 ergs cm" 2 s" 1 ) 

5.9 ±0.5 

6.1 ±0.6 

AF (10" 17 ergs cm" 2 s" 1 ) 

120 ± 50 

120 ± 60 

Lbright (10 29 ergs s" 1 ) 

24.2 ±0.2 

24.5 ±0.4 

Lfaint (10 29 ergs s" 1 ) 

1.15 ±0.09 

1.2 ±0.1 


Note. — Maximum and minimum luminosities of C0M-M71A as derived from the model light curve. 
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Fig. 1.— Postfit timing residuals for M71A, as functions of orbital phase and date. The 
dashed line indicates orbital phase 0.25, the PSR superior conjunction. The asymmetric 
eclipse spans just over 15% o f the orbit wh ich is typical for BW systems at radio frequencies 


below ~ 2 GHz (e.g., Fruchter e t al. 1 9881 ). 
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Fig. 2.— Variation of the time of passage through ascending node (orbital phase 0), com¬ 
puted for overlapping segments of data and holding all other timing parameters fixed at their 
nominal values. 
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Phase 


Fig. 3.— An observation of M71A on MJD 52798 (2003 June 8), using 3 WAPPs to cover 
300 MHz of contiguous bandwidth. The abrupt disappearance of the PSR at the start of 
eclipse, as well as the slight dispersive delay on reappearance, are clearly visible. The cumu¬ 
lative pulse profile is plotted twice at the top of the figure. 





27 



Fig. 4.— Primary dataset HST images of the 5" x 5" region around the nominal position 
of M71A. The filters and the orbital phases are labeled in each panel. The solid circle is 
centered on the radio position and it has a radius of 0.2" (which is larger than the formal 
uncertainty from PSR timing). The dashed circle is centered on the X-ray counterpart and 
it has a radius of 0.5". The relatively bright star on the left border of the solid circle is 
the candidate optical companion proposed bv iHuang et al.l ( 20101 ). COM-M71A is clearly 
visible inside the solid circle in the right panels (corresponding to the inferior conjunction 
of the PSR, where the companion reaches maximum brightness), while in the left panels (at 
superior conjunction of the PSR) it is below the detection threshold. 
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0.5 1 1.5 

Orbital Phase 


Fig. 5.— Light curves of COM-M71A in the F606W and F814W Liters separately (upper 
and middle panels) and for the combination of the two (bottom panel), obtained after a 0.95 
mag shift of the F814W magnitudes. All curves are folded with the radio parameters and two 
periods are shown for clarity. Circles mark the observed points, arrows are the magnitude 
upper-limits for the images where the star is below the detection threshold. The black curve 
in each panel is the best analytical model obtained from the combined light curve and then 
adapted to each Liter. 





29 



Orbital Phase 


Fig. 6.— Light curves of the candidate companion proposed by iHuang et ah (2010), folded 
with radio orbital parameters. The absence of any magnitude modulation as a function of 
the orbital phase is the definitive confirmation that this object is not connected to M71A. 
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Fig. 7.— CMD of M71 with highlighted in blue the region occupied by COM-M71A during 
the whole orbital period, as predicted by the light curve model (see text, Figure [5] and 
Table 3). 
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McOM (M 0 ) 


0.2 0.08 0.05 0.04 0.035 0.0325 



Fig. 8.— Reprocessing efficiency of the PSR emitted energy as a function of the inclination 
angle, assuming three different values of the Roche Lobe filling factor and a PSR mass of 
1.4M 0 . The thickness of each strip corresponds to the range of spin-down energies measured 
for this PSR (see text). The horiz o ntal dashed line at r/ = 15% is a typical reprocessing 
efficiency reported in Breton et al. ( 2013 ). On the top axis, the companion masses for a 
PSR mass of 1.4M 0 are reported. 
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0.5 1 1.5 

Orbital Phase 


Fig. 9.— Optical light curves of COM-M71A (this work; black points and lines) and COM- 
M5C (from Pallanca et al. 2014a,, gray points and dashed lines), with magnitudes reported 


to absolute values. The gray solid lines are COM-M71A model adapted to COM-M5C to 
reproduce the observed points. 
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Fig. 10.— CMD of M71 (black dots) and M5 (gray dots). The blue shaded region is the 
position of COM-M71A along the whole orbital phase, as derived from the light curve model 
(see Figure]?]). The red po int and the gray area are the indicative position of the COM-M5C 
(see Pallanca et al. 2014a, for more details). Both the objects are located between the MS 
and the WD cooling sequence, suggesting common properties of these two BW companions. 











